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the quark content of baryons and mesons. The level of this review is intended for undergraduates or 
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through deep inelastic scattering; structure functions and what they reveal about proton structure; 
and prospects for further studies with new and upgraded facilities, particularly a proposed electron- 
ion collider. 
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I. INTRODUCTION 

These lectures were given at Aligarh Muslim University, Aligarh India, at a nuclear physics workshop attended 
by Indian graduate students and postdocs. These lectures provide a pedagogical review of the past 50 years of 
experimental and theoretical studies of the quark and gluon structure of the nucleon. Since many speakers at this 
workshop will discuss structure functions extracted from deep inelastic scattering reactions, we define these quantities 
and review what can be extracted from them. We hope that specialists in this field will tolerate our brief and 
superficial summary of many years' study of strongly-interacting systems. 




FIG. 1: Richard H. Dalitz, 1925-2006. 

These lectures are dedicated to my mentor, Richard Dalitz, who was from 1963-2000 Royal Society Professor of 
Physics at Oxford University and who made seminal contributions in particle phenomenology. Dalitz pairs (electron- 
positron pairs from decay of a high-energy photon) and the Dalitz plot were named after him. In addition he 
made important contributions to our understanding of strange particles and hypernuclei, and to particle spectroscopy 
through the constituent-quark model. I appreciated his insights into theoretical physics, his keen intuition about 
physical problems, and his encouragement of young physicists. 

II. DISCOVERY OF NEW PARTICLES WITH DEDICATED ACCELERATORS 

The development of both accelerators and detectors progressed rapidly through the first half of the 20 th century 
PQ. Ernest Lawrence's cyclotron ushered in a new era which provided both much higher beam energies and more 
intense beams than were previously possible. This was followed by the development of synchrotrons and mastery of 
magnetic focusing techniques, which allowed particles to be contained in a circular pipe surrounded by individual 
focusing magnets. Shortly following World War II, development began on a succession of machines based upon the 
synchrotron principle. Advances in accelerator technology were accompanied by corresponding advances in particle 
detectors. The bubble chamber was an extremely important tool in discovering new particles. In particular strange 
particles like the S and f2 were discovered in bubble chambers. Another major detector advance was the development 
of the spark chamber. 
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Most of the particles discovered before the 1950s were either completely stable (so far as we know), like the electron 
or proton, or had relatively long lifetimes. For example, the lifetimes of charged pions and kaons are roughly 1CP 8 
seconds, and the lifetimes of the A and charged £ are of order 10 -10 seconds. When charged (neutral) particles move 
with sufficiently large velocities, they leave measurable tracks (gaps) in emulsions or bubble chambers between their 
production and decay vertices. 

Starting in the early 1950s, researchers began to discover states which were much more short-lived than these earlier 
particles. These states were extracted from scattering experiments at what were then the highest available energies. 
Fig. [2] shows elastic and total cross sections for ir + — p scattering. Note the very prominent peak at CM energy 
just near 1.2 GeV. This is what was then called the 'P33 resonance' in ir — p scattering. The notation for this state 
denoted the angular momentum (the spectroscopic notation P for L = 1), and the isospin / and total spin J in the 
notation L^i^j- In current terminology this is the A(1232) resonance [2]. Since the total width for this state is of 
order 100 MeV then from the Heisenberg uncertainty principle which can be written as A_EAt > h, we can infer that 
the lifetime of the A(1232) is of order 10~ 23 seconds. 
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FIG. 2: Total and elastic X-sections for tv + — p scattering, as a function of the laboratory beam momentum in GeV/c (note 
logarithmic scale), also the CM energy in GeV. 

For such short-lived states the production and decay vertices could never be separated, regardless of the particle's 
velocity. As scattering energies increased, and as nuclear targets were bombarded with different probes, more and more 
peaks appeared in cross sections. Determining the properties of these resonances became a major area of research 
beginning about 1951, and spectroscopic studies of resonance production have continued ever since. Information about 
these resonances, and comparative ratings of the reliability of claimed resonant states, are compiled in the biennial 
report of the Particle Data Group [2]. 



A. Resonances and Scattering Amplitudes 

Here is an elementary review of bound states and scattering resonances [3]- Consider a two-body scattering process; 
for simplicity we will initially assume spinless non-relativistic particles in an S- wave state, interacting via a potential 
V(r). For the time being we will consider two-body scattering in only the elastic channel. The radial wave function 
can be written as R(r) = ip(r)/r, and the radial Schrodinger equation satisfied by i]j{r) for a given value of the relative 
momentum k is 

-^^ + U{r)Mr) = k 2 Ur) (1) 

In Eq. 0, U(r) = 2mV(r), m is the reduced mass of the two-body pair, and we adopt units where %=1. 

Strong interaction potentials are short-ranged, so we expect U(r) — > as r — > 00. In this case at large r we can 
write 

lim Mr) = a(k)e- lkr + b(k)e lkr (2) 
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We can consider the wavefunction ip to be an analytic function of k. A sufficiently strong attractive potential may 
support one or more bound states. These are confined states that decay exponentially. They will occur as discrete 
solutions of the Schrodinger equation of Eq. ([1]) corresponding to imaginary momentum k = iks- For this particular 
value of the momentum, the coefficient a{iks) in Eq. ([2]) must vanish; in this case the wavefunction will decay 
exponentially at large distances, i.e. V>fc B ( r ) ~ e~ knT ' . 

The properties of bound states and resonances are related to the behavior of the scattering amplitude in the complex 
k plane, or equivalently the complex E plane. The non-relativistic relation between energy and momentum is 



k 2 /(2m) + E th 



(3) 



From Eq. (J3J) it is clear that as the phase of the complex variable k varies from — > 27r, the phase of E will change 
from — > Air. Thus one sheet of the complex k plane will map onto two sheets of the complex E plane. For the time 
being consider only the elastic channel (inelastic channels will produce additional sheets of the complex E plane) . 

A bound state corresponds to a pole of the scattering amplitude T{E) in the complex k plane at positive imaginary 
value iks- From Eq. (|3]), bound state poles will appear on the real E axis, below the continuum threshold Eth- The 
upper sheet of the complex E plane is called the physical sheet (the shaded sheet on the left in Fig. [3]), and the lower 
sheet is an unphysical sheet. The two sheets are joined by a cut running along the positive real E axis, beginning at 
continuum threshold. 
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FIG. 3: Left: the complex E plane, the physical sheet (top, shaded) and an unphysical sheet (bottom). A cut runs along the 
real E axis, beginning at continuum threshold. A resonance is shown as an 'X' on the unphysical sheet. Right: X-section vs. 
energy for a Breit-Wigner resonance, whose amplitude is given by Eq. Q. 



If one crosses the cut starting from above, moving from the physical sheet onto the unphysical sheet, the T matrix 
is continuous. However, T is discontinuous if one crosses the cut beginning and remaining on the physical sheet. 
The discontinuity of the T matrix across the cut plays an important role in constructing dispersion relations for the 
scattering amplitude [3]. 

Imagine that one has an attractive potential strong enough to produce a bound state. As one decreases the strength 
of the attractive potential, at some point it will no longer support a bound state. The bound state pole may move 
onto the second sheet of the complex E plane and become a resonance. A resonance generally produces a peak in the 
scattering cross section at some energy above elastic threshold. The prototypical non-relativistic resonance amplitude 
was posited by Breit and Wigner [3] , who wrote the scattering amplitude in the vicinity of the resonance in the form 



T(E) = g(E)[E- M + iT/2]- 



(4) 



This leads to a peak in the cross section vs. energy for E ~ M, a so-called 'Breit-Wigner' resonance. The right-hand 
side of Fig. [3] plots the qualitative features of the resulting cross section, which is proportional to the absolute square 
of the scattering amplitude. 
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Eq. Q shows that the resulting cross section will have a peak at energy E — M; furthermore, if the quantity g(E) 
is relatively constant in the vicinity of the peak, then the full width at half maximum around the peak will be given 
by the quantity T. So the location and width of the resonance peak are characterized by the two quantities M and T. 
Eq. Q shows that the scattering amplitude will have a pole in the complex plane at E = M — iT/2. The responance 
pole lies on the unphysical sheet of the complex E plane, below the real E axis, as indicated on the left in Fig. [3| 

We can also examine the properties of resonances by defining the scattering or S- matrix |3J. The T matrix is 
defined as the amplitude of an outgoing wave solution, corresponding to a unit incoming wave. If we restrict ourself 
to the elastic channel, the relation between the S and T matrix is 

T(E) = S( - E } - 1 ; S(E) = e 2lS ^; T{E) = e l5{E) sin(S(E)) (5) 
2i 

Unitarity requires that the ^-matrix have absolute value one. An attractive (repulsive) potential tends to produce 
positive (negative) phase shifts. In a one-channel situation with a Breit-Wigner resonance, Eq. @ shows that T (E) 
will be purely imaginary at the resonance position E = M. By inspection of Eq. ([5]), a purely imaginary T-matrix 
corresponds to 8{E = M) = ir/2. The location of an elastic resonance (with no background) is characterized by a 
phase shift of 90°. 
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FIG. 4: Argand diagram. At low energies the T matrix moves on the unitary circle. The phase 5 = it/2 at an elastic resonance. 
When inelastic channels open, the T matrix moves inside the unitary circle. 

An Argand diagram represents a convenient graphical way to study the behavior of the scattering amplitude, 
through a two-dimensional plot of the T matrix vs. energy. The horizontal axis is the real part of the T matrix and 
the vertical axis is the imaginary part. For purely elastic scattering the T matrix moves on a circle (the unitary 
circle) of radius 1/2 centered at the point (0, 1/2). The origin of the Argand diagram corresponds to 6 — 0. The angle 
subtended by any point on the unitary circle (relative to the origin) is 26. This is illustrated in Fig. [4] which shows 
an Argand diagram for n — r\ scattering [5]. For sufficiently low energies, only the elastic channel is open, and the T 
matrix moves on the unitary circle. The phase shift 6 — n/2 at the position of the ao(980) resonance. The smaller 
the width T, the faster the phase shift will move through the resonance. 

At higher energies, inelastic channels open up. The complex E plane must be expanded to include inelastic channels. 
Every time an inelastic channel opens, a new cut appears on the real E axis at the inelastic channel threshold, and 
new unphysical sheets appear. Locating a scattering resonance becomes more difficult, as one must determine both 
the resonance position and the sheet on which it occurs. In the presence of inelastic channels the scattering amplitude 
becomes a matrix. The quantity Tfi(E) is defined as the coefficient of an outgoing wave in channel / corresponding 
to unit incoming wave in channel i. Unitarity requires that in the presence of inelastic channels, the clastic T matrix 
element must lie inside the unitary circle, and Fig. [^demonstrates this. Finding resonances in multichannel reactions 
is likely to rely on sophisticated fitting of scattering amplitudes, or on models that allow one to analytically continue 
the T matrix in the complex E plane. 

We should emphasize that the simple Breit-Wigner resonance is almost never observed in actual experiments. The 
7r — 7T resonance p{77Q) is one of the few examples where this simple resonance picture is valid [2 . In most other cases 
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the situation is far more complex. The goal in analyzing resonant reactions is to identify the location of the pole in 
the complex energy plane, as this pole is the fundamental dynamical quantity. 

B. Identification of Resonances with the Dalitz Plot 

In the 1950s the catalog of resonances grew at first modestly, then exploded. From a few resonances, one rapidly 
accumulated a dozen, then eventually greater than one hundred resonant states! A major contributor to this 'zoo' 
of resonant states lay in the ability for accelerators to produce secondary beams of pions and kaons. These could be 
scattered from protons and nuclei, and the final states could be analyzed. 

Our simple analysis of resonances was carried out for spinless particles. For particles with spin it is necessary to 
combine spin and orbital angular momentum in order to determine the total spin and parity of a given resonance. 
Sophisticated multi-channel partial-wave analyses were developed to extract the scattering amplitudes from experi- 
mental reaction data. As energies increased many final states involved three or more strongly interacting particles. A 
major advance in resonance phenomenology was the introduction of the Dalitz plot, a method to analyze three-body 
(and higher) final states from scattering reactions [5J. 




FIG. 5: Dalitz plot for the reaction p + p — > 37r . The shading represents the frequency of events, with blue denoting few and 
red signifying many events. Four resonances are visible in this plot. 

To understand the Dalitz plot, consider a reaction a + b — > 1 + 2 + 3 leading to three final-state particles. For the 
Dalitz plot one graphs the data vs. the squared invariant mass of one pair of final-state particles on one axis and the 
invariant mass of a second pair on a perpendicular axis. The data points fall inside well defined kinematic limits, 
since one can show m\ 2 + wif 3 + ^23 = const. Each observed event corresponds to a point inside the Dalitz plot. If 
the scattering state decays uniformly into three particles, the resulting points will be uniformly distributed. However, 
two-body resonances appear as enhancements in the Dalitz plot, at the location of the resonant mass. Furthermore, 
the structure of these bands depends on the spin and parity of the resonant state (s). 

Fig. [5] shows the Dalitz plot for the reaction p + p — > 3ir°, from the Crystal Barrel experiment [7]. The data points 
are represented by shading, with blue denoting few and red signifying many events. The axes are chosen as the 
invariant mass of two different ir° pairs. The pion Bose symmetry requires that the Dalitz plot be symmetric about 
any axis that exchanges two pions. Four resonances are identified in this plot. The spin-0 /o(1500) appears as a 
simple enhancement. Two spin-2 resonances, the /2(1270) and /2(1565), have a more complicated shape. The final 
resonance, the /o(980) with spin 0, also appears on the plot. Because of quantum mechanical interference between 
amplitudes, the /o(980) appears as a dip rather than a peak in the distribution of events. The Dalitz plot has proved 
incredibly useful in enabling experimentalists to locate resonances in reactions with multi-particle final states. In 



G 



addition, this technique provides a quick estimate of the spin and parity of the observed state. 

III. THE QUARK STRUCTURE OF HADRONS 

By the beginning of the 1960s, there were of order 100 known resonances. These could be classified into baryon 
resonances, strongly interacting states with half-integral total angular momentum J, and meson resonances, states 
with integral (or zero) spin. What had initially been exciting discoveries of a few new states had become almost an 
embarrassment. This apparent over-abundance of resonances was for awhile quite perplexing. The situation seemed 
to cry out for introduction of a classification scheme that would make sense of the masses and spins of the plethora 
of observed resonances. 

The first significant breakthrough came in 1960-61, when Nc'cman and Gell-Mann independently showed that nearly 
all mesons and baryons could be grouped in multiplets defined by the symmetry SU(3) [8]. In this scheme, dubbed 
the 'Eightfold Way' by Gell-Mann, the observed resonances occurred in families when plotted vs. their electric charge, 
isotopic spin and strangeness values. The graph at left in Fig. [6] shows the SU(3) classification picture for the baryon 
octet (a multiplet of spin-parity J p = 1/2+ particles) and decuplet with J p = 3/2 + . The right-hand figure shows 
the extension of octet and decuplet multiplets when one includes charm quarks [2j . 
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FIG. 6: Left: the baryon octet and decuplet as described by SU(3) symmetry. The vertical axis is strangeness, the horizontal 
axis is the third component of isotopic spin, and lines of constant electric charge are indicated. Right: extensions of the octet 
and decuplet including charm quarks. 

The baryon octet consists of a pair of particles (the neutron and proton) with zero strangeness, four particles with 
S = — 1 consisting of an isotropic triplet (the Es) and an isosinglet A, and a doublet of S with strangeness S = —2. 
One could derive formulae relating the masses of particles within a multiplet; these formulae were generally accurate 
to within a few percent. Significant support for this classification scheme came in 1962 with the discovery of the S* 
doublet with strangeness — 2 This left only one 'missing' member of the baryon decuplet, and Gell-Mann predicted 
a new baryon with spin-parity 3/2 + and strangeness —3 [8] dubbed the multiplet mass formulae predicted a mass 
1680 MeV. In 1964 the was discovered at Brookhaven in bubble chamber photographs of K~ + p reactions [10] ; 
its mass was estimated as 1686 MeV. 

Similar SU(3) multiplets, with corresponding mass formulae for the multiplet members, could be constructed for 
the known pseudoscalar and vector mesons. The next step was to find an underlying theory that explained why 
mesons and baryons should inhabit such multiplets. This answer was provided in 1964 independently by Gell-Mann 
|llj and Zweig |12j . who proposed that baryon and meson structure could be explained naturally by the existence of 
quarks (Zweig called them aces). The idea was that there existed three types of quarks, now called flavors (up, down 
and strange), each with its corresponding antiquark. Baryons were composed of three quarks and a meson was a 
quark-antiquark pair. Individual quarks were fermions with intrinsic spin 1/2; strangeness was an additive property, 
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with the up and down quarks having S = and the strange quark S = — 1. However, in order to produce the observed 
spectrum of particles, it was necessary that quarks possess electric charges that were a fraction of the proton's charge, 
specifically the up quark had charge +2e p /3 and the down and strange quarks had charge — e p /3 (an antiquark had 
opposite strangeness, charge and intrinsic parity from its respective quark). 

A. Simple Quark Models of Hadron Structure 

Following the introduction of quarks, the proton was considered to be a three-quark combination consisting of two 
up quarks and one down quark. This would produce the requisite electric charge of +1 and strangeness zero. In a 
similar vein, the K + meson with strangeness S = +1 and charge +1 could be composed of an up quark and a strange 
antiquark. Furthermore, the spin and parity of a baryon or meson could in principle be constructed by forming various 
spin-parity combinations of the requisite quarks. 

Several groups began to consider what are now known as constituent quark models to investigate the apparent 
properties of these quark combinations, in terms of the known spins, parities and masses of known states |13j . 
Beginning in the mid 1960s, such efforts had some notable successes. Many of the properties of the baryon octet and 
decuplet could be reproduced in a constituent quark model with up and down quarks of roughly 340 MeV mass and 
a strange quark mass in the vicinity of 450 MeV. The three valence quarks in the proton would occupy the lowest 
states in a confining potential such as an harmonic oscillator. The up quark spins were parallel to the direction of 
the proton's spin, and the down quark was anti-parallel, with the three quarks coupled to overall spin 1/2 and isospin 
1/2. The A ++ (1232) state was composed of three up quarks, in the same confining potential. Here the quark spins 
are aligned, producing overall spin-parity 3/2 + . The strength of the color magnetic interaction could be inferred from 
the N — A mass splitting. The A is composed of uds valence quarks. The N — A splitting gives the mass of the 
strange constituent quark. In this model of the A, the light quarks are coupled to spin and the spin of the strange 
quark is parallel to the A spin direction. In a similar vein, the lowest nonets of pseudoscalar and vector mesons are 
made by coupling qq pairs. 

Despite their apparent success in describing the lightest hadron multiplcts, constituent quark models encountered 
a series of problems. Some of these were: (1) The constituent quark model predicts that the first excited baryons 
should have negative parity. However, the lowest-lying nucleon excited state, the iV(1440), has positive parity [2"1[TB]. 
Furthermore, the masses, spins and parities of excited baryonic and mesonic levels disagreed with constituent quark 
predictions. (2) Constituent quark masses consistent with the lowest baryon and meson multiplets dramatically over- 
predicted the pion mass of about 140 MeV. (3) The constituent quark model suggested that quarks in a nucleon were 
rather lightly bound and non-relativistic. Arguing from experience with atomic nuclei, the constituent quark model 
suggested that it should be reasonably easy to eject an individual quark from a proton and observe it in a detector. 
(4) This suggested that fractional charges could be emitted from nucleons. Many extremely sensitive experiments 
were mounted to search for fractional charges; apart from one or two anomalous results, no isolated fractional charges 
were ever detected [14] . (5) Constituent quark models required the quarks to be totally symmetric upon interchange 
of space, spin and isospin quantum numbers. Constituent quark (fermion) wave functions thus violated the Pauli 
principle. 

B. Development of Quantum Chromodynamics 

The non-observation of free quarks and violation of the Pauli principle led some to postulate that quarks might 
not be real but might simply represent a mathematical 'mnemonic' device enabling one to predict meson and baryon 
properties. However, several arguments against the 'reality' of quarks were soon overcome. First was the discovery 
of the property of color, an additive property of quarks which takes on three values [15]. All stable particles exist 
only in combinations of quarks and antiquarks that couple to zero net color. For baryons, the resulting three-quark 
wavefunction is completely color-antisymmetric. Thus quark wavcfunctions obey the Pauli principle when one includes 
color degrees of freedom. 

The next significant evidence for quarks came from deep inelastic scattering (DIS) reactions [16] . These involved 
scattering of leptons (charged leptons, or neutrino-induced charge-changing reactions) from nucleons at very large 
momentum transfers; one observed only the outgoing lepton and no information regarding the struck nucleon(s). 
Such reactions showed distinctive features consistent with scattering from elementary pointlike objects inside the 
proton and neutron. Later in these lectures, we will discuss in considerable detail the history of DIS reactions, what 
they reveal about baryon structure, and our current understanding of quark-parton distributions in the nucleon. 

The quark model was extended with the discovery of three additional quark flavors. In 1974 the charm quark was 
discovered, with mass roughly 1500 MeV [17] . In 1977 the bottom quark was discovered, with a mass about 4500 
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MeV [18 . Finally 1995 marked the discovery of the top quark, with a mass roughly 171,000 MeV [15] . The quarks 
can be grouped into three doublets (u,d) (c,s) and (t,b); in each pair the first quark has charge +2/3 and the second 
quark —1/3. Except for the light quarks, the doublets are also characterized by exceptionally large mass splittings. 

With the discovery of heavier quark flavors one expects to see additional multiplets of particles containing the heavy 
quarks (the right side of Fig. [6] shows low-lying multiplets with 4 flavors) . We now know that the up and down quark 
masses are much smaller than the 'constituent' quarks (u and d 'current quark' masses are roughly 4 MeV and 9 MeV, 
respectively). Thus, the nucleon mass contains very large contributions from sea quarks and gluons, and there is no 
way that a non-relativistic perturbative treatment can explain the properties of baryons containing only light quarks. 
On the other hand, charm and bottom quarks are sufficiently massive that simple models should accurately describe 
hadrons containing these quarks. Heavy quark effective theory (HQET) has made great use of the simplifications that 
arise due to the very large mass of the heavy quarks (20] . 

We have now a detailed theory of quark interactions. This has culminated in a non-Abelian gauge theory of 
quarks interacting through the exchange of colored gluons, the carriers of the strong force between quarks. This 
theory is known as quantum chromodynamics or QCD [21 . A first breakthrough was an understanding of the 
phenomenon of asymptotic freedom. Through the Nobel Prize- winning work of Gross and Wilczek and Politzer |22j , 
it was demonstrated that QCD had the property that at high energies (or alternatively, at very short distances) 
the interaction between quarks became progressively weaker and weaker. At these very high energies, one could 
approximate interactions with quarks in terms of free quark interactions. Asymptotic freedom demonstrated why 
deeply inelastic scattering at very high momentum transfers looks like scattering from free constituents. 
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FIG. 7: Left: the value of the QCD strong coupling constant a s (E) vs. the energy E in GeV. The decrease of the coupling 
constant with energy demonstrates the phenomenon of asymptotic freedom. Right: schematic picture of quark confinement. 
As quarks are separated, a 'string' or flux tube forms. When the string is broken a new q — q pair is formed. 

The left-hand graph in Fig. [7] shows the strong coupling constant a s as a function of energy. The property of 
asymptotic freedom is demonstrated by the monotonic decrease in a s with increasing energy. Conversely, the strong 
coupling strength increases rapidly at low energies. This means that an understanding of QCD at low energies must 
deal with very strongly-interacting systems. This is additionally complicated by the realization that gluons have 
strong self-coupling, unlike the situation in QED. QCD offered formidable challenges to the development of models 
of strongly interacting systems. 

A second major feature of QCD is that an individual colored quark experiences a force that grows monotonically at 
large distances. This is shown schematically in the right-hand picture in Fig. [7] As a quark-antiquark pair is separated 
a string or flux tube forms. If one attempts to break the string confining an individual quark, eventually one will 
produce an additional q — q pair (a meson), rather than liberating the quark. Confinement has never been proved 
analytically in QCD, however it provides a natural explanation for the experimental failure to observe individual 
quarks. Confinement also occurs in lattice gauge QCD. 

The QCD Lagrangian has the form 

C QC o = q^V^-m)q- 1 -G^ 
G^u = d^G® — dyG^ — gf abc G b ^G^ (6) 
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In Eq. (|6|) one sums over all flavors of quark fields, represented by q, the quantities G° represent the eight gluon fields, 
and theierm is the QCD analog of the electromagnetic tensor F^ v which occurs in quantum electrodynamics. 
An excellent summary of QCD is given in the Handbook of Perturbative QCD [23] , 

As we have mentioned, the problem of solving the QCD Lagrangian at low energies turned out to be formidable. 
Our best estimate of u and d quark mass is less than 10 MeV, which leads to the astonishing conclusion that the 
'current quark' masses make up only about 2% of the proton mass! The remainder of the mass must consist of very 
large contributions from gluons and the sea. At present, the only viable means to solve QCD at low energies appears 
to be lattice gauge techniques. In this case one discretizes space-time on a lattice. Good numerical results can be 
achieved by rotating to Euclidean (imaginary) time. The continuum limit is achieved as the lattice spacing a — > 0. 
One discretizes the Yang-Mills action for the QCD Lagrangian, defines the quark fields on the lattice sites, and uses 
Monte Carlo techniques for solving the resulting multi-dimensional integrals. Lectures on lattice gauge theory will be 
given at this workshop by N. Mathur [21]. 



IV. DIS AND THE PARTONIC STRUCTURE OF BARYONS 



Our current view is that hadrons are composed of structureless elementary particles, quarks and antiquarks, whose 
strong interactions are mediated by exchange of colored gluons. We now review the history of the experimental data 
that demonstrated the quark nature of baryon structure, and indicates that quarks have no internal structure. The 
first significant evidence for the elementary nature of quarks came from deep inelastic scattering (DIS) reactions |16j . 
These involved scattering of leptons (initially charged leptons, or neutrino-induced charge-changing reactions) from 
nucleons at very large momentum transfers. In these reactions one observed only the outgoing charged lepton and no 
information regarding the struck nucleon(s). Such reactions carried out initially with electron beams at the Stanford 
Linear Accelerator Center (SLAC) and with neutrinos at CERN, showed distinctive features consistent with scattering 
from elementary pointlike objects inside the proton and neutron. Initially these constituents were called 'partons' but 
these were later shown to be quarks and antiquarks |13j . 

Consider DIS reactions induced by scattering of charged leptons from a nucleon, mediated by exchange of a virtual 
photon, as shown schematically in Fig. [8] We can describe the process in terms of two relativistic invariants, which 
we will define as Q 2 and x, defined as Qr = —q 2 and x = Q 2 / '(2Mu); here the 4-momentum transfer Q 2 is large, and 
v = E — E' is the lepton energy loss. 




FIG. 8: Schematic diagram of deep inelastic scattering of charged lepton, with initial and final momenta k and k' , from a 
nucleon with momentum P, through exchange of a virtual photon. 



hadron tensor W^, where 



The cross section for DIS with a charged lepton can be written as the contraction of a lepton tensor and a 



dxdy Q 4 ^ 
W, v = ^- g ^ + q ^F^ x ,Q 2 ) + ^F^x,Q 2 ) (7) 

The cross section for lepton-nucleon DIS (mediated by photon exchange) thus depends upon two structure functions 
and Each structure function depends on the invariants x and Q 2 . In a frame of reference where the proton is 
moving with infinite speed, the quantity x gives the fraction of the nucleon's momentum carried by the struck parton. 
In Fig. [8j one could also exchange a weak vector boson Z between the charged lepton and nucleon. This would give 
rise to a third structure function denoted F3, which would be associated with parity- violating transitions. 
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Over the past thirty years, many experiments have been devoted to measuring these structure functions. The left- 
hand side of Fig. [9] shows measurements of the structure function F^. Values of F^ are plotted for various values of x, 
vs. Q 2 . The curves for different x values are offset so that they do not lie on top of one another. The values have been 
obtained from a series of different experiments [25] . There are electron-induced experiments from the SLAC linear 
collider and the HERA electron-proton collider (its two detectors, HI and ZEUS), and experiments using high-energy 
muons from the BCDMS and New Muon Collaboration (NMC) at CERN and experiment E665 from Fermilab. 
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FIG. 9: Left: the structure function F% obtained from DIS from charged leptons on protons. F2 is plotted for different values 
of a; as a function of Q 2 . Right: polarized structure function xgi(x,Q 2 ), obtained for a fixed value Q 2 — 10 GeV 2 and plotted 
for p, D and n. 

The first dramatic feature of the resulting structure function is the precision that has been obtained with these 
F2 measurements. Except for the very highest values of Q that can be achieved for a given x, the errors are very 
small indeed. The second striking feature is that, over a Q 2 range of five orders of magnitude, the structure function 
F% varies very little. The constancy of the Fi structure functions was clear evidence that scattering was occurring 
from elementary pointlike constituents; otherwise the structure function would be governed by a form factor, and over 
this wide range of Q 2 a form factor should decrease by orders of magnitude. Friedman, Kendall and Taylor shared 
the 1990 Nobel Prize in physics for their DIS measurements in e — p reactions at SLAC that first demonstrated this 
'scaling' behavior in the structure functions [25] , 

One can also measure polarized structure functions. For scattering of unpolarized leptons on polarized nucleons, 
the cross section asymmetry resulting from virtual photon exchange is related to the polarized structure function 
gi(x,Q 2 ), where one has 

1 =*' ~ x 9l (x, Q 2 ) (8) 

ax ay 

In Eq. ([8]), the cross section a(N^) is taken for a nucleon whose spin is parallel to the direction of the incident lepton 
beam, and subtracted from the cross section cr(N^) where the nucleon's spin is antiparallel. The right-hand side 
of Fig. [9] plots xgi vs. x for experiments with protons, deuteron, and neutron (the neutron result is extracted from 
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experiments on 3 He) [37]. Here, the experiments are all extrapolated to the same value Q 2 = 10 GeV 2 . One can 
also extract additional structure functions when the nucleon spin is transverse to the plane of the four-momentum 
transfer. 



A. A Qualitative Picture of Proton Structure 



In Sec. |III A| we presented a naive picture of the structure of the proton, based on the constituent quark model. 
Based on our current understanding of QCD dynamics, and our knowledge that the 'bare' masses of u and d quarks 
are extremely small, we need to modify this simplistic picture. We know that at low energies quarks are coupled to 
gluons with a coupling constant. Furthermore, gluons exhibit self-coupling and are also coupled to quark-antiquark 
pairs. These couplings are shown schematically in the left-hand figure in Fig. |10| 




FIG. 10: Left: schematic picture of quark and gluon couplings. A quark can radiate a gluon; gluons exhibit self-coupling; and 
a gluon can radiate a quark-antiquark pair. Right: a cartoon showing a 'valence plus sea and glue' picture of the proton. 



Imagine that we start with a nucleon containing just three 'valence' quarks. Through the QCD couplings shown 



in the left-hand side of Fig. 10 the valence quarks can radiate gluons, which in turn can produce more gluons and 
quark-antiquark pairs. Starting with a proton that contains only 'valence' quarks, QCD radiation can produce a 
proton containing a 'sea' of gluons and quark-antiquark pairs, as depicted by the cartoon on the right in Fig. [10] 

We have presented two naive pictures of the proton, one with just three constituent quarks and the second with 
many additional sea quarks and gluons. What we expect in a given DIS reaction will depend upon the value of 
x. Note that from the left-hand side of Fig. [10] quark or gluon radiation involves sharing parton momentum with 
the additional constituents. As QCD radiation produces more and more partons, each parton will carry less of the 
proton momentum and will correspond to partons with progressively smaller values of x. Consequently, the quantity 
x which denotes the fraction of the proton's momentum carried by the struck quark, will determine which picture of 
the nucleoli will be 'seen' in that reaction. At large values of x the nucleon will look like three 'valence' quarks; as x 
decreases we expect a 'sea' of quark-antiquark pairs and gluons. 

A crucial ingredient is the probability of finding partons that carry a given fraction x of the proton's momen- 
tum. Quantities that satisfy this definition are called Parton Distribution Functions, or PDFs. Define the quantity 

q(x), q = u,u,d,d, For a proton with very high momentum, q(x) is the probability that a quark of that flavor 

carries a momentum fraction between x and x + dx; similarly a gluon PDF g(x) gives the probability that a gluon 
carries momentum fraction x. We should also include a dependence upon Q 2 , as we will later show that PDFs have 
a slow (logarithmic) dependence on Q 2 . 

From Fig. |10| antiquarks in the proton will arise mainly from gluon radiation of q — q pairs. Thus we expect all 
antiquarks to be part of the nucleon 'sea' that is radiated by gluons. Light (u and d) quarks have two sources: one 
component of these PDFs is part of the 'valence' distribution; a second component arises from gluon radiation. Now 
q and q arising from gluon radiation have equal probability. For a given quark flavor we thus define valence quark 
distributions q v (x) = q(x) — q{x). If all antiquark distributions arise from gluon radiation, then subtracting the 
antiquark PDF from the quark PDF cancels out the 'sea' quark contribution leaving just the 'valence' part. 



B. Relation Between DIS and Quark PDFs 



In a preceding section we showed that DIS reactions could be described in terms of a small number of structure 
functions. An obvious question is: what is the relation between the structure functions and the quark PDFs that we 
have just defined? To lowest order in QCD, structure functions have very simple relations in terms of PDFs. For 
example, the unpolarized and polarized structure functions F% and gj that occur in DIS from charged leptons (arising 
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from virtual photon exchange between electron and proton as shown in Fig. [8]) have the form 



F](x, Q 2 ) = e) [q 3 {x) + qj{x)] = 



4(u + u) + d + d+s + s 
9 



g J(x,Q 2 ) = £]^[A<fc(aO+A&(a 



4(Aw + Am) + Ad + Ad + As + As 
18 



(9) 



In Eq. (|9j|, ej is the charge of a quark with flavor qj, and Aq(x) = ~q(x) — *q~(x), where ~q(x) (*q(x)) is the probability 
for finding a quark whose spin is parallel (antiparallel) to the proton spin. Eq. ^ includes contributions from three 
flavors. 

If we include DIS from charged-current v leptoproduction, e.g. e~ + p — > v + X, then the reaction is characterized 
by a virtual W~ being absorbed by the proton. In that case three of the structure functions that appear in this 
reaction can be written in lowest order QCD as 



xF™-(x,Q 2 ) 



2x [u + d + s] 
2x [u — d — s] 
Au + Ad + As 



(10) 



The additional structure function xF^ is present because exchange of the weak vector boson W~ contains a piece 
that violates parity ; he re the structure function F% characterizes the parity-violating part of this interaction. Note in 
both Eqs. ^ and (10), we have not included the dependence of the PDFs on x and Q 2 . 

The above equations show that different observables are sensitive to different combinations of PDFs. Thus by carry- 
ing out a scries experiments employing lcpton, neutrino and hadron beams, and including unpolarized measurements 
with both longitudinal and transversely polarized beams and targets, one can gain differential sensitivity to specific 
quark and gluon PDFs. By amassing a large body of experimental data from many different experiments one can 
perform a global analysis of all of the experimental data and extract best values for the quark and gluon PDFs [2] . 

We mentioned previously that PDFs at different Q 2 are related by QCD, through what are called the DGLAP 
evolution equations [SB]. These are convolution relations relating the PDFs for a given value of x at different Q 2 
values. Thus if one knows the PDFs at one (sufficiently large) value of Q 2 , they can be determined at higher values of 
Q 2 through these evolution equations. The non-singlet (valence) PDFs satisfy their own integro-differential equation, 
while the sea quark and gluon distributions are coupled, 
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In Eq. (Ill, the quantity nj refers to the number of active quark flavors participating in the reaction, and the evolution 



p. 



can be determined from perturbation theory 



V. EXPERIMENTALLY DETERMINED PARTON DISTRIBUTION FUNCTIONS 

Over the past three decades a very large number of precision experiments have been performed, designed to elucidate 
the partonic structure of nucleons. The pioneering experiments at SLAC were fixed-target experiments, with a beam 
of high-energy electrons incident on a fixed target. Such experiments can reach the largest values of x, and are being 
continued today at Jefferson Laboratory. Another class of experiments involved colliding beams. For example, the 
HERA accelerator at DESY, which has recently been closed down, was an asymmetric collider where 30 GeV electrons 
collided with 820 GeV protons. Fig. [9] shows that the two HERA detectors HI and ZEUS covered an enormous range 
of x and Q 2 , and also reached to the lowest values of x. HERA also had a smaller fixed-target program, HERMES, 
which studied electron collisions with polarized nucleoli targets. 
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Useful information regarding parton distributions has been obtained from Drell-Yan processes [29] ■ These processes 
involve the production of large invariant-mass /i + — \i~ pairs from nucleon-nucleon and nucleon-nucleus reactions. 
These processes occur when a q (q) from the projectile annihilates a q (q) of the same flavor from the target, producing 
a virtual photon that decays to a pair of oppositely charged muons. Such processes are quite sensitive to the q PDFs 
in the proton. 

Information on s quark distributions is obtained from neutrino charge-changing reactions. In such processes an 
incident produces a \x~ and a virtual W + . The W strikes an s quark and produces a c quark. The charm 
quark subsequently undergoes semi-leptonic decay producing a /i + . The signature of these reactions is production of 
oppositely charged high-energy muons from a neutrino beam, and provides the most reliable information on s quark 
PDFs. Production of /i + — fi~ pairs from a v beam is sensitive to s PDFs [3"U] . 

Information from a wide array of high-energy experiments is placed in a global fit to parton distribution functions. 
From these experiments one can extract the structure functions. By now, the relation between the structure functions 
and PDFs is presented at least to next to leading order (NLO) in QCD, and sophisticated global fitting routines 
are employed to extract the parton distributions. Two of the best known of these are employed by the CTEQ 
collaboration [31], and by the MRST group [32] • The DGLAP evolution equations of Eq. (Ill are used to connect 
parton distributions at different Q 2 values. The resulting PDFs are plotted at one particular value of Q 2 ; PDFs at 
different values of Q 2 can be obtained from the evolution equations. 
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FIG. 11: Proton PDFs. Values of xq(x) determined from global fits to high-energy data evolved to Q 2 = 10 GeV 2 . Left: 
width of band gives a measure of the uncertainty in the PDF. Yellow: determined from HI HERA data; green: determined 
from ZEUS data. Right: a linear plot in x, separating different sea quark flavors. 

Fig. [TT] shows our current knowledge of quark PDFs. Quark and gluon PDFS multiplied by x are plotted vs. x. 
The left figure is a semi-log plot and the right-hand figure is a linear plot that separates different sea quark flavors. 
The PDFs are evolved to Q 2 — 10 GeV 2 . The HERA data is most extensive. There are small differences between 
the PDFs depending on data from the HI detector (yellow band) or from ZEUS (green band) [25]. The width of the 
bands gives an indication of the uncertainty in the PDFs. The quantity xS is the sum of all sea quark PDFs. 

Our naive picture suggested that at large x, the proton should be composed primarily of valence quarks. This agrees 
with the analyses; for large x the valence quarks dominate the PDFs. For small x the sea and glue distributions rise 
rapidly. Note that in the left graph of Fig. [TT] the sea and gluon distributions have been divided by 20. Thus the sea 
and glue grow extremely rapidly at small x. 

The valence PDFs obey quark normalization conditions. Their first moments give the total number of valence 
quarks in the proton, (u v { x)) = 2, (d v (x)) — 1. where the brackets denote integration over all x. One thus expects 



that u v (x) <~ 2d v (x). Fig. 11 shows that this is approximately the case. 



In the right-hand side of Fig. 11 the sea PDFs arc separated by flavor. The s quark distribution is roughly half the 



average of u and d sea quark PDFs. The c quark distribution i s m uch smaller than the other sea quark distributions, 



due to the much larger c quark mass. The right side of Fig. 11 shows that the d distribution is larger than the u 
distribution. We can understand this from 'meson cloud' models of the nucleon |33j . Such models include contributions 
particularly from the pion cloud. The proton can undergo virtual transitions p — > n + n + . The valence quark content 
of the 7r + is u — d. Thus one contribution to the d distribution arises from scattering from the d in the 7r + cloud of 
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the proton. 'Meson cloud' models can account at least qualitatively for the magnitude and x distribution of the d — u 
distribution in the proton. 

VI. SURPRISING FEATURES OF PARTON DISTRIBUTION FUNCTIONS 

In the preceding section we reviewed qualitative features of quark PDFs obtained from global fits to high energy. 
Now we will focus on some of the surprises that have occurred as we have unraveled the structure of the proton. 

A. Momentum carried by the Proton 

In our simplest picture with a proton composed of three constituent quarks, each 'constituent' quark should carry 



roughly 1/3 of the proton's momentum. The valence quark distributions shown in Fig. 11 peak at lower values of 
x; u v peaks at x ~ 0.2 and d v peaks at an even smaller x value. A major surprise occurred with the quantitative 
understanding of the distribution of the proton's momentum. The total proton momentum can be represented by a 
sum rule 

l^xiq^x) + qj{x)) + xg{x)^ =1 (12) 

Eq. ( |12[ ) expresses the fact that the proton momentum can be divided into the total amount carried respectively by 
quarks, antiquarks and gluons. At high values of Q 2 , the total proton momentum carried by valence quarks is roughly 
35%, sea quarks carry 15%, and gluons carry 50% of the momentum. This differs dramatically from the naive picture 
where valence quarks would carry all of the proton's momentum. It is one more indication of the extremely important 
role of glue, both in constituting the mass of the nucleon and now in terms of the proton linear momentum. 

B. Nuclear Modification of Structure Functions 

Another major surprise is the realization that even at very high energies, structure functions in nuclei differ from 
those in the proton. First let us discuss what we expect for parton distributions in the neutron. The operation of 
charge symmetry (CS) interchanges p and n; this is a specific rotation in isospin space (a rotation of 180° about the 
'2' axis in isospace) that interchanges p and n labels. At low energies, charge symmetry is obeyed extremely well, 
with most low-energy amplitudes obeying CS at the 1% level or better. So, it is 'natural' for CS also to hold at high 
energies. In principle one can test the validity of parton charge symmetry; no violations have been observed and 
present upper limits of partonic charge symmetry are at the level of several percent |34j . 

If one assumes the validity of charge symmetry at the partonic level, then all PDFs for the neutron can be written 
in terms of those for the proton. Charge symmetry predicts that u n (x) — d p (x) and d n (x) = u p (x); it also predicts 
that s and c PDFs should be identical for n and p. Analogous equations hold for antiquarks. Since heavier nuclei have 
roughly equal numbers of neutrons and protons, in comparing nuclear structure functions to those in the nucleon it 
makes sense to compare heavier nuclei with the deuteron. 

To date, there have been extensive measurements of the F 2 structure function for DIS of charged leptons from 
nuclei. One plots the F2 structure function per nucleon arising from virtual photon exchange, and constructs the ratio 
of the structure function for a nucleus with A nucleons with that for the deuteron, 

R A (x) = F 2 A (x)/F 2 D (x) (13) 



Fig. 12 plots the quantity Ra(x) vs. x for various nuclei, from Arneodo [35] . The top curve shows the qualitative 
behavior of Ra{x) vs. x. The bottom points are the results of this ratio for experiments on four nuclei. The diamonds 
and open circles are results on C and Ca respectively from NMC. The solid squares show results on Al from SLAC. 
The solid triangles are results on Xe from E665. 

The experimental points show rather dramatically that the respective structure functions F 2 have a distinct A 
dependence. Furthermore, the nuclear effects on the F 2 structure functions can be divided into roughly three distinct 
regions: a 'shadowing' region for small x < 0.05; a 'Fermi smearing' region for very large x; and what is called the 
'EMC effect' in the region of x between the shadowing and Fermi smearing regions. Fig. 12 shows dramatically that 
even at high energies, there remain significant differences between DIS from nuclei and scattering from individual free 
nucleons. 



15 




FIG. 12: The ratio of Fi structure functions for nuclei to those for the deuteron, vs. x. Top: a curve that shows the qualitative 
behavior of this ratio vs. x. Bottom: results of this ratio for experiments on four nuclei. 



The ^4-dependence of Ra{%) for intermediate x is called the 'EMC effect' after the definitive experiments in the early 
1980s from the European Muon Collaboration at CERN |36j . These experimental results were extremely surprising. 
At least in the regime of intermediate x, it was widely assumed that at sufficiently high energies, nuclear DIS would 
look just like DIS from quarks in individual free nucleons. 

Since x = Q 2 /(2Mu), the very large x region corresponds to the largest values of Q 2 for a given incident energy. 
On a nuclear target, scattering of a virtual photon at large x values will be greatly affected by the nuclear Fermi 
momentum. For collisions at very large Q 2 to the nucleus, nucleons with Fermi momentum moving towards (away 
from) the virtual photon will experience collisions at a relatively much larger (smaller) value of Q 2 . Since F2 — > 
rapidly at large x, Fermi smearing effects will produce nuclear F2 that are much larger than that for a free nucleon. 
Thus we expect to see a very rapid increase in Ra(x) at sufficiently high x due to these Fermi smearing effects, in 
agreement with the results shown in Fig. [T2| 



1. Nuclear Effects in the Shadowing Regime 



Fig. 12 shows that in the region of small x < 0.05, the ratio Ra{x) decreases monotonically with decreasing x, 
and Ra{x) shows some sign of approaching a constant value at very small x. The experimental points show that 
this curve is not universal but has a characteristic A dependence, with the 'shadowing' corrections becoming larger 
with increasing A. For a given incident energy the small x region corresponds to the smallest values of Q 2 . Even 
for very high energies, sufficiently small x values can correspond to the range 1 < Q 2 < 5 GeV 2 . In this region, 
viewed from the rest frame of the nucleus, we expect to sec effects arising from transitions of the virtual photon to 
vector mesons. From the CVC (Conserved Vector Current) hypothesis, a virtual photon will occasionally make a 
transition to a vector meson, particularly a p(770). Since photons interact weakly with matter we expect the impulse 
approximation to be valid, in which case Ra{x) should be one. However, the p interacts very strongly and at these 
energies will be strongly absorbed. To a good approximation the cross section for the p should vary as the area of 
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the nucleus, a p ~ R 2 ~ ^4 2 / 3 . These are called 'shadowing' corrections. As soon as a p encounters a nucleon it is 
absorbed. Nucleons on the 'far' side of the nucleus are 'shadowed' from seeing the p by the nuclcons on the 'near' 
side. 

The nuclear F 2 structure function per nucleon arising from 7* — > p transitions should behave like A 2 / 3 /A ~ A^ 1 / 3 . 
One would expect to see Ra(x) ~ A -1 / 3 if 100% of the nuclear cross section arose from such transitions. Although 
quantitative calculations are significantly more complicated than this simple picture, the qualitative behavior is correct. 
First, shadowing corrections decrease F 2 in heavier nuclei and one expects Ra(x) < 1; the shadowing corrections 
are largest at the smallest values of x; and finally shadowing corrections become larger with increasing A. All of 
these are observed in the data. If we consider nuclear DIS reactions in the infinite momentum frame, then shadowing 
should occur through gluon recombination and interference effects. A picture due to Gribov [57] relates shadowing to 
diffractive effects. This picture relates two different scattering phenomena and leads to interesting predictions but it 
is not yet clear that it provides quantitative agreement with experiment. 



2. The 'EMC Effect' 

The final nuclear effect involves a slight increase Ra(x) > 1 just below x = 0.1, followed by a monotonic decrease 
in the region 0.1 < x < 0.6; this is called the 'EMC Effect' from its discovery in the early 1980s by the European 
Muon Collaboration. Until that time one expected that nuclear DIS in this region would just be the sum of individual 
lepton-nucleon DIS; thus we would expect F 2 (x) per nucleon to be the same regardless of A. 

The EMC Effect appears to be nearly universal; the curves of Ra(x) for various nuclei lie almost on top of one 
another. The EMC Effect inspired more than 1000 papers that considered the origin and explanation of this nuclear 
effect [36]. Kulagin and Petti [38] claim that three different mechanisms play a role in the EMC Effect. The first two 
involve nuclear binding and off- shell effects on nuclear structure functions. Treating bound nucleons non-relativistically 
and using weak-binding approximations, one can write F 2 as a product of free F 2 structure functions for p and n, 
convoluted with the virtuality of the bound nucleons. An additional contribution to F 2 comes from scattering of the 
virtual photon from the meson cloud, particularly the pion cloud. The left side of Fig. 13 plots the ratio Ra(x) vs. x 
for gold [35]. The dotted curve includes only Fermi momentum plus nuclear binding (FMB) effects. The dashed curve 
includes the FMB plus off-shell (OS) effects. The dot-dashed curve adds nuclear pion (PI) effects; and the solid curve 
is the full result including nuclear shadowing (NS) effects. Shadowing is significant only in the shadowing regime. 
Nuclear pion effects are significant only for x < 0.2. In the EMC region the most important processes are claimed to 
be binding and off-shell effects, although others have argued that nuclear pion effects are dominant [3"6] . 



C. Polarized DIS and the Spin of the Proton 



As was mentioned in Sec. IV B one can extract the quark longitudinal spin-dependent PDFs xAu(x) and xAd(x) 
through measurements using polarized protons. Eq. ([8| relates the polarized structure function gi(x) to the asymmetry 
in deep inelastic scattering for unpolarized leptons on longitudinally polarized nucleons. Eq. ^ relates the polarized 
structure function to the spin-dependent PDFs. On the right in Fig. [13] we plot the longitudinal spin-dependent PDFs 
that have been extracted from experiment. 

The experimental results [39 are from the Spin Muon Collaboration (SMC) (open squares) and HERMES (solid 
circles). HERMES took semi-inclusive DIS measurements from polarized proton or deuteron targets, where one 
measured a final-state meson (ir or K) in coincidence with the scattered lepton. Such semi-inclusive DIS experiments 
allow one to enhance contributions from particular spin-dependent quark flavors. As can be seen from Fig. |13| one 
obtains Au v (x) > and Ad v (x) < 0, and in magnitude Au v (x) > |Ad v (a;)|, in agreement with naive quark model 
expectation. It is also clear that the sea quark contribution to the nucleon spin is small and consistent with zero for 
all measured x values. 



1. The 'Proton Spin Crisis' 



We can write a sum rule for the total spin of the proton, 

1 _ 1 

2 ~ 2 



AT, + AG + L q + L g , (14) 



In Eq. (14 1, AS is the total spin carried by quarks and antiquarks, AG is the total spin carried by gluons, and the 
quantities L„ and L g are respectively the total orbital angular momentum carried by quarks and gluons. The capital 
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FIG. 13: Left: the ratio R A {x 
of Ra{x)\ bottom: semi-log plot 



Curves evolved to Q 2 = 2.5 GeV 2 . 



from Eq. ( 13 1 of F2 structure functions for gold to those for the deuteron 
Right: spin-dependent PDFs. Top: xAu v (x) vs. x; middle 



vs. x. Top: linear plot 
xAd v (x); bottom: xAq 3ea (x). 



letters refer to the first moment over x, i.e. 

AS = (J2 Ag(x) + Aq(x)) ; AG = (Ag(x)) (15) 

In the constituent quark picture the first moments of quark spin PDFs are given by 

AU = -, AD = -i AS* = 0, AS = AU + AD + AS = 1; (16) 



Eq. (16 1 shows that in the naive quark picture all the proton spin arises from the quarks, with zero contribution from 
gluons or orbital angular momentum. 

The first quantitative measurements of the quark spin PDFs gave results like AS ~ 0.2 ± 0.2. Not only was the 
total proton spin carried by quarks far smaller than 100%, these results were consistent with zero proton spin being 
carried by quarks! This surprising and puzzling result was termed the 'proton spin crisis.' We have now obtained 
more precise results and the current best value is AS ~ 0.3 ± 0.05. So the total proton spin carried by quarks is 
somewhat larger than it had appeared, and it is no longer consistent with zero. Nevertheless, one very important 
quantitative question is exactly how do we account for the spin carried by the proton [40 ? A related question is 
whether or not the current experimental result AS ~ 0.3 for the proton spin carried by quarks is surprising. 



2. Contribution of Polarized Glue to Proton Spin 

Previously we saw that gluons carried a surprisingly large fraction of the proton momentum. It seems natural to 
ask whether gluons might also carry much of the proton spin. In order to answer this important question, dedicated 
experiments to measure gluon contributions to the proton spin have been mounted at three laboratories. There is 
the COMPASS experiment at CERN, experiments from the RHIC Spin group at the Relativistic Heavy Ion Col- 
lider (RHIC) at Brookhaven, and the HERMES experiment at HERA. We will discuss particularly the COMPASS 
experiment. 

At COMPASS, the goal was to measure gluon contributions to the proton spin through photon-gluon fusion |41j . 
In this process, shown schematically in Fig. |14[ a polarized muon beam couples to a virtual photon. This in turn is 
coupled to the polarized gluon distribution in the proton through production of a quark-antiquark pair. The final 
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scattered muon is observed in coincidence with a final hadronic state which has strong overlap with the q — q pair. 
The COMPASS experiment focused on two different hadronic states; the first was 'open charm,' observation of a final- 
state hadron containing a charm quark. The second situation involved two final-state hadrons with large transverse 
momentum p T . 




FIG. 14: Left: schematic diagram illustrating 'photon-gluon fusion.' The virtual photon from a polarized muon couples to 
the gluon distribution through production of a quark-antiquark pair. Right: Ag(x) / g(x) vs. x. Star: COMPASS open-charm 
measurement; circle and downward-facing triangle: COMPASS high-p T measurements with different Q 2 cuts; square: SMC 
high-p r measurement; upward-facing triangle and diamond: two different HERMES frigh-p T measurements. 



The right-hand side of Fig. 14 shows measurements of Ag(x) / g(x) from several laboratories |42j . The star represents 
the COMPASS open-charm measurement. The circle and downward-facing triangle are COMPASS high-p T measure- 
ments corresponding to different Q 2 cuts. The square is a high-p T measurement from the Spin Muon Collaboration 
(SMC) group. The upward-facing triangle and the diamond are two different HERMES high-p T measurements. The 
data is all for values 0.05 < x < 0.2. The values of Ag(x)/g(x) are all small and in most cases the errors in both Ag 
and the average value of x are relatively large. The errors are sufficiently large that the sign of Ag is not clear. The 
two curves show fits to these data with AG > and AG < 0. 

If one takes the COMPASS values one obtains for the total gluon contribution to the proton spin, AG/G = 
-0.57±0.41 (stat) ±0.17 (syst) for the open charm measurements, and AG/G = 0.016±0.058 (stat) ± 0.0.055 (syst) 
for the high-p T measurements [43]. Combining the COMPASS measurements with the SMC, HERMES and RHIC 
spin results, the best value of AG is small; however, it is still statistically possible that AG could be large. Fig. [14] 
does not show RHIC spin results that give comparable values for AG/G in a similar x region. 



3. Current Status of the Proton Spin Question 

At the present time, despite the large errors on the gluon contribution to the proton's spin, it appears that polarized 
gluons are unlikely to carry the majority of the proton's spin. We are left with the question of just what fraction of 
the spin of the proton is carried by polarized glue, or by quark or gluon orbital momentum. An area that is being 
pursued at the moment is hard exclusive processes, such as deeply virtual Compton scattering or DVCS (where a 
charged lepton interacting with a nucleon exchanges a virtual photon with high Q 2 , producing a real photon in the 
final state), or exclusive meson production in lepton-nucleon scattering. These processes can be described in terms of 
so-called generalized parton distributions or GPDs. It was shown by Ji that a particular moment of the GPDs could 
be related to the total orbital angular momentum of the quarks |44j . 

Investigations of DVCS have been undertaken at HERA and at Jefferson Laboratory. A major focus of the 12 
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GeV upgrade at Jefferson Lab will be studies of DVCS and exclusive meson production. The physics of the 12 GeV 
upgrade will be reviewed in detail in the first lecture by Tony Thomas at this workshop [15] . Thomas' second lecture 
will summarize his recent work with Myrher on the spin of the proton, which will address the questions raised in this 
section. 



VII. NEW AND PROPOSED FACILITIES TO STUDY NUCLEON STRUCTURE 

To date, much information on the partonic structure of the nucleon has been obtained from a series of high-energy 
machines. The first experiments that showed scaling in high-energy deep inelastic scattering were carried out at 
SLAC. Following these measurements were continuing fixed-target experiments at SLAC, and both fixed-target and 
collider studies at CERN and Fermilab. The HERA asymmetric electron-proton collider proved to be a nearly ideal 
facility to map out parton distributions over a very wide range of x and Q 2 . Neutrino and antineutrino beams at 
Fermilab and CERN allowed one to separate out contributions from different quark flavors. Studies of Drell-Yan 
processes at Fermilab and CERN provide processes sensitive to sea quark distributions in the nucleon. 

Now that HERA has ceased operations, experiments at CERN and Fermilab continue. Various new facilities or 
upgrades of existing ones have the capability of extending our knowledge regarding the partonic structure of the 
nucleon. The first is the Large Hadron Collider or LHC which will very shortly begin operation at CERN. The 
primary thrust of LHC operations will be first to find the Higgs particle, a new massive particle which is predicted to 
appear as the result of the spontaneous breaking of electroweak symmetry. Next, the LHC will focus on discovering 
cither supersymmetric particles or other new phenomena beyond the Standard Model. However, the LHC also could 
have potential applications in determining the partonic content of hadrons. 

There are also several new facilities, or upgrades of existing facilities, which could examine the partonic content of 
hadrons. Here I will discuss one upgrade of an existing facility and one proposed new facility in this regard. The first, 
the upgrade of the electron accelerator at Thomas Jefferson Laboratory, is the subject of the first lecture by Tony 
Thomas at this workshop [45] . As a result, I will provide only a brief summary of this project; I refer the reader to 
his talk on this subject. The second is a proposed new electron-ion collider. Two possible versions of this machine 
have been suggested, one which would be located at Brookhaven National Laboratory and a second which would be 
the subject of a future upgrade at Jefferson Laboratory. 

I will very briefly discuss the 12 GeV upgrade of the CEBAF accelerator at Thomas Jefferson National Laboratory 
|46j . The upgrade will provide 12 GeV continuous electron beams to a new Hall D experimental area, which will 
produce real photon beams that should provide precision spectroscopy. This should help to study possible 'exotic' 
states, particularly excited mesons that contain one or more gluons. In addition, the upgrade would provide up to 11 
GeV beams to the existing halls A, B and C. 

As was mentioned in the preceding section, a major focus of effort following the Jefferson Lab upgrade would 
be on hard exclusive processes such as DVCS or exclusive meson production. It has been shown that these hard 
exclusive processes can be described in terms of integrals over quark distributions. These are termed 'generalized 
parton distributions' or GPDs. GPDs can provide information regarding the longitudinal momentum and transverse 
position of quarks in the nucleon, and could test the Ji sum rule [44 . Following the upgrade, Jefferson Lab will have 
unique kinematic capabilities. It will be able to access parton distributions in regions of x and Q 2 that complement 
those accessible at laboratories like CERN or HERA. In particular, after the upgrade Jefferson Lab will be able to 
reach regions of Bjorken x ~ 0.6 and values Q 2 ~ 8 GeV 2 . 

A proposed new facility is an electron-ion collider. This would consist of an electron beam colliding with either 
light ions or heavy ions. Plans to date call for achieving a value s = 20 — 100 GeV 2 , with high luminosity. There 
are two suggestions for such a facility |47j . In each case this would involve building a new accelerator to complement 
an existing machine. The first case would involve building an electron ring at Brookhaven, to collide with hadrons 
from the current RHIC facility. The left-hand side of Fig. [15] shows a schematic picture of such a facility. A proposed 
electron ring would be added to the RHIC complex at Brookhaven; the ring would provide electrons with energies in 
the range 5-10 GeV, with the possibility of a recirculating linac injector. The electrons would collide with beams of 
heavy or light ions from the RHIC accelerator. 

The second scenario would involve building a light ion accelerator at Jefferson National Laboratory. This would 
produce hadrons of 50-100 GeV, to collide with electrons of roughly 5 GeV from CEBAF. The right-hand side of Fig. [15] 
shows a schematic picture of such a facility. A hadron ring (in this diagram, shown in a figure-eight configuration) 
would be added to the CEBAF electron accelerator at Jefferson Lab. The ring would provide protons and light 
ions with energies in the range 50-100 GeV. The proposed hadron accelerator would feature four intersection points. 
Present estimates suggest that very high luminosities (£ ~ 10 34 /cm 2 -sec) might be achieved in such a collider. 

With an electron-ion collider one could provide quantitative measurements of quark spins through semi-inclusive 
DIS. It would have the possibility of measuring PDFs (and particularly gluon distributions) at very small Bjorken 
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FIG. 15: Conceptions of possible electron-ion colliders. Left: eRHIC. An electron ring added to the RHIC complex at BNL. 
Right: electron-light ion collider at Jefferson Lab. A hadron ring would be added to the CEBAF accelerator, accelerating 
protons and light ions. 



x. Gluon probabilities in the nucleoli increase extremely rapidly at very small x values, as can be seen from Fig. |11| 
The large nucleoli densities in the center of the nucleus could in principle support extremely large gluon densities. 
If gluon distributions continue to increase in the interior of the nucleus, at some point nuclear gluon probabilities 
would exceed the unitarity bound. Thus one expects at some point for many-body effects to slow down this increase. 
Some have suggested that this will lead to a qualitatively new cooperative effect termed a 'color-glass condensate' 
[48j . An electron-ion collider, particularly eRHIC which supports heavy-ion hadronic beams, would be able to probe 
this dynamical region. 

An electron-ion collider would also have the capabilities of exploring generalized parton distributions over a fairly 
wide kinematic region. Such colliders could also provide detailed quantitative studies of the structure of mesons, 
particularly pions. One other area that could be studied with an electron-ion collider is the process of fragmentation. 
This is the process whereby a quark produces a final hadron. Possible kinematics for electron-ion colliders would have 
the capability of exploring both beam fragmentation and the so-called target fragmentation regions. 

The author was supported by NSF contract PHY-0555232, and he acknowledges useful discussions with A.W. 
Thomas and A. P. Szczepaniak. 
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